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without any free parameter. The effective running coupling is small in the IR and the approximation 
can be improved by inclusion of higher order terms. 
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It is well known that perturbation theory (PT) breaks 
down in the infrared (IR) limit of QCD. The coupling 
grows up and eventually diverges at the low-energy 
scale Aqcd ~ 200 MeV, so that our knowledge of 
the very important low-energy phenomenology, includ¬ 
ing hadron masses and quark confinement, must rely on 
non-perturbative methods and mainly on numerical lat¬ 
tice simulations. In the IR, we still miss a full analytical 
description that could be regarded as sound and reliable 
as PT is in the ultraviolet (UV). On the lattice, by non- 
perturbative renormalization schemes like moment sub¬ 
traction (MOM), the renormalized coupling turns out to 
be finite and not too large in the IR, so that its divergence 
must be regarded as an artifact of standard PT. A finite 
coupling would suggest that a perturbative approxima¬ 
tion could be set up and that the failure of standard PT 
is a consequence of the usual perturbative approach more 
than a problem with the general perturbative method by 
itself. 

In this paper, by an unconventional setting of the per¬ 
turbative method, PT is shown to be working in the IR 
limit of pure SU(3) Yang-Mills (YM) in Landau gauge, 
yielding a first-principle analytical description that repro¬ 
duces the lattice data fairly well, without any free param¬ 
eter that was not in the original Lagrangian. Moreover, 
the approximation can be improved, order by order, by 
inclusion of higher order terms, as for any other pertur¬ 
bative expansion. 

A key role is played by the dynamical mass that the 
gluon acquires in non-perturbative approximations, as its 
appearence is prohibited by gauge invariance at any order 
of standard PT. That is a signature of an intrinsic failure 
of the usual perturbative expansion in the IR and the 
gluon mass can be regarded as a pure non-perturbative 
effect. The success of any perturbative expansion re¬ 
lies on the choice of the zeroth order approximation that 
should capture most of the physics, leaving minor details 
to be defined by inclusion of the interaction. While a 
free massless gluon is a good starting point in the UV, 
no finite mass can arise in PT if it is not already present 
in the zeroth order of the expansion. On the other hand, 
it has been recently shown that, by inclusion of a mass 


term in the Lagrangian, PT turns out to be perfectly 
viable in the IR, yielding one-loop results that fit the lat¬ 
tice data very well[I|, [3|. Thus the gluon mass seems to 
absorb most of the non-perturbative effects and having 
included the mass in the zeroth order propagator, the ef¬ 
fects of the interaction can be treated by standard PT. 
Moreover, in the UV limit the mass becomes irrelevant, 
and the modified Lagrangian maintains the correct high 
energy limit. 

However, inclusion of a mass by hand would not be 
a satisfactory solution of the problem for many reasons. 
First of all, the mass term is forbidden by gauge invari¬ 
ance, together with any counterterm that should be re¬ 
quired for its renormalization. In addition, while a free 
mass-parameter would be useful for interpolating exper¬ 
imental and lattice data, it would be an extra parameter 
that was not in the original Lagrangian, yielding a phe¬ 
nomenological model at most. 

The choice of a more clever zeroth order Lagrangian 
does not require the inclusion of extra terms because we 
have the freedom of splitting the original action in any 
two parts, yielding different perturbative expansions. In 
fact, that freedom suggests a way to optimize the expan¬ 
sion by a variational evaluation of the best zeroth order 
action[3-l6|. In other words, we may add a mass term to 
the zeroth order Lagrangian and subtract the same term 
in the interaction, without affecting the total Lagrangian 
that maintains its gauge invariance. While the total ac¬ 
tion does not depend on that mass, any finite-order ex¬ 
pansion would depend on the mass parameter that can be 
optimized by some variational ansatz in order to get the 
best perturbative expansion. The optimized PT is not 
gauge invariant at any finite order, but gauge invariance 
must be recovered approximately in the physical observ¬ 
ables if they converge, order by order, towards the exact 
ones. 

For a pure SU(N) YM theory in Landau gauge, drop¬ 
ping color indices, the zeroth-order gluon propagator can 
be taken as 
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Figure 1: First and second order graphs contributing to the 
ghost self energy and the gluon proper polarization. The filled 
circles are the bare vertices of Lint while the cross is the mass 
counter-term m^. 


where Aq is the propagator of a massive free gluon 

Ao(p)”^ =+ TO^ (2) 

while the zeroth-order ghost propagator Gq remains un¬ 
changed (massless). By standard Dyson equations, the 
dressed gluon and ghost propagators read 

Mp)~^ = -p^ + - n*(p) (3) 


G(p)-i=p2-E*(p) (4) 


where If* is the proper polarization and E* is the proper 
self energy. These functions are to be evaluated, order 
by order, by PT as expansions in powers of the total in¬ 
teraction, including a mass counter-term that cancels the 
mass in the total Lagrangian. Namely, the total interac¬ 
tion Ltot can be written as 


^tot - ^ 




( 5 ) 


where 5 is a power-counting parameter to be set to one 
at the end of the calculation and Lint is the original in¬ 
teraction of YM theory, including the ghost vertex and 
three- and four-gluon vertices that contain powers of the 
bare YM coupling g. Thus, the perturbative expansion 
is not a loop-wise expansion in powers of the coupling 
g, but must be written as a 5-expansion where graphs 
with a different number of loops may coexist at the same 
orderQ. In Fig.l the one-particle irredicible graphs are 
displayed up to second order in S. 

Of course the proper functions S*, H* do contain di¬ 
vergences and a regularization scheme is in order. While 
dimensional regularization is suited for standard loop- 
wise expansions where gauge invariance is preserved at 
any finite order, here gauge invariance is already lost and 
a more physical non-perturbative scheme would be pre¬ 
ferred, like a simple cutoff. Restricting all integrals to 
p 2 ^ ^2 Euclidean space, the graphs in Fig.l can 

be evaluated for any choice of the bare YM coupling g 
and for any value of the parameter . Moreover we can 
make contact with lattice simulations where a natural 


cutoff is provided by the lattice spacing and we can use 
the same non-perturbative multiplicative renormalization 
scheme, thus avoiding any insertion of counter-terms in 
the Lagrangian. As for a lattice, we assume that the 
bare coupling g = g{A) is a function of the cutoff and the 
physical content of the theory is invariant for a change 
of scale A ^ A' followed by a change of coupling g ^ g' ■ 
Then, physical renormalized functions that do not de¬ 
pend on the cutoff nor on the coupling can be defined at 
an arbitrary scale p by the usual multiplicative renormal¬ 
ization. That is equivalent to say that all bare functions, 
at different couplings, should go one on top of the other 
by scaling, yielding a test for the accuracy of the pertur¬ 
bative approximation. 

A drawback of the simple choice of a cutoff is the occur¬ 
rence of spurious quadratic divergences besides the usual 
logarithmic terms that are observed on the lattice. These 
quadratic terms must cancel in the exact propagators and 
polarization functions so that their appearance can be re¬ 
garded as an artifact of the finite-order approximation. 
Moreover, while the exact functions cannot depend on 
m?, the finite-order expansion does depend on the choice 
of the mass parameter and the spurious divergences are 
sensitive to the choice of m?. That would suggest a sim¬ 
ple criterion for determining the best mass parameter: 
we assume that the best expansion is the one where the 
quadratic divercences are vanishing (or smaller). At any 
finite order, that choice would give an expansion free of 
spurious divergences and no free parameters at all, allow¬ 
ing for a challanging comparison of its predictions with 
lattice data. Besides its formal elegance, that choice is 
quite resonable if one wants to recover the correct scaling 
properties of the theory. 

An other similar choice could rely on the use of Steven¬ 
son’s principle of minimal sensitivity^ that would re¬ 
quire to make the physical observables stationary. In 
fact, while in principle vn? can be chosen at will, any sen¬ 
sible physical prediction should not depend on its choice. 
The mass parameter is not a physical observable by 
itself, and different strategies for its optimization may 
lead to very similar physical results if m? is close to a 
stationary point where small changes in the mass of the 
zeroth order propagator are compensated by the changes 
on the interaction. 

It is not difficult to show that the second-order proper 
functions in Fig.l can be made free of quadratic diver¬ 
gences by a special choice of the mass parameter. We 
may neglect the last graph of Fig.l, namely the two-loop 
sunrise graph, as it has been shown to be very small and 
basically constant Q. Its inclusion would just make the 
derivation awkward without any substantial change in 
the result. Moreover, we may neglect subleading terms 
that vanish in the limit A ^ oo. The quadratic diver¬ 
gences can be extracted by the polarization function at 
p = 0 as the difference n*(p) —n*(0) contains logarithmic 
terms at most. By an explicit calculation of the graphs 
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in Fig.l we 


findQ 
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where the effective bare coupling a = 9A^a/(167r) and 
a = (;^/(47r). For TV = 3 the effective coupling a « 
is quite small even when g is rather large. We recognize 
a term (a^ log^) arising from the two-loop tadpole: when 
A ^ m the infinite series of tadpoles would eventually 
require a resummation of the large logarithms. 

The quadratic divergence cancel exactly if the mass 
parameter is taken as 



exp 




( 7 ) 


which has the evident non-perturbative behaviour ^ 
exp(—const./a) in the limit a —>■ 0. Inserting Eq.(I7]) back 
in Eq.® the term (a^ log^) cancels exactly, yielding the 
simple linear expression 


n*(o) = 


. 5 16 , 
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By insertion of Eq. ([7]) again, we can make the logarithmic 
divergence implicit through the coupling and write 


n*(p) = w? 


' 64 1 ' 


+ [n*(f.) - n*(0)| (9) 


where the last term 11* (p) — 11* (0) is just given by the 
sum of gluon and ghost one-loop sunrise graphs, the only 
ones that depend on p. By Eq. ([3|) the dressed propagator 
acquires the dynamical mass n*(p)] that in the limit 
p —^ 0 defines a physical mass 

+ 1 °) 

which survives after the cancellation of the mass terms in 
the zeroth order propagator and in the interaction. After 
renormalization, the mass is expected to be invariant 
for any change of coupling and cutoff. 

In order to compare with the phenomenological model 
of Ref. [3| we explore the optimized second-order expan¬ 
sion for = 3 in the large-coupling range g = 4—5 where 
that model was found in perfect agreement with the lat¬ 
tice data. We find a good agreement with that model, 
which is not surpising as Eq.([9]) would be equivalent to 
the result of a one-loop calculation with a phenomeno¬ 
logical mass inserted by hand and set equal to M^. 

Eirst of all, in Fig.2 we check that the gluon propaga¬ 
tors can be renormalized by a multiplicative scale factor 
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Figure 2: Log-Log plot of the renormalized gluon propagator 
for A = 3 in arbitrary units for g — 4, 4.4 and 4.9. 


and translated one on top of the other by a translation 
in a log-log plot. As one would expect by a finite-order 
approximation, the scaling properties are not stisfied ex¬ 
actly but we can still extract a renormalized propagator 
which is almost independent of the coupling. The resid¬ 
ual dependence suggests that the agreement with lattice 
data would be at its best for a special value of the cou¬ 
pling that turns out to he g ^ 4.9 in the phenomeno¬ 
logical model of Ref.[^. That small dependence should 
decrease by inclusion of higher order terms. 

Having no experimental data on pure SU(3) we fix the 
physical ener^ scale by a comparison with the lattice 
data of Ref. [9] and renormalize the propagator at the 
same energy of that work. The cutoff turns out to be 
A = 1.05 Gev at g = 4.9. That value is small but not 
too far from the inverse lattice spacing that is 1/i ^ 



Figure 3: The renormalized gluon propagator in physical units 
(A = 1.05 Gev) for A = 3 and g = 4.9. The points are the 
lattice data of Ref.[^ for P = 5.7 and L = 96. 
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Figure 4: The renormalized ghost dressing function F{p) = 
p^G{p) for N = 3 and g = 4.9 in the same physical units 
of Fig.3 (A = 1.05 Gev). The points are the lattice data of 
Ref.0 for l3 = 5.7 and L = 96. 

1.16 GeV in the simulation^. The renormalized gluon 
propagator is shown in Fig.3 together with the lattice 
data. The agreement is very good for a first-principle 
second order calculation with no free parameters. The 
physical mass from Ea. dTOll turns out to be M = 0.56 
GeV and is constant in the range g = A — 4.9 as shown in 
table I. That value is not far from the phenomenological 
parameter m = 0.54 GeV extracted by a best fit of data 
in Ref.j^l- 

In Fig.4 the renormalized ghost dressing function 
F{p) = p^G{p) is reported with the same energy units 
(A = 1.05 Gev) and compared with the lattice data. The 
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Figure 5: The running coupling as defined in Eg. (Hill , for 
N = 3 and g = 4.9, with the same renormalization constants 
of Fig.3 and Fig.4 (solid line). The points are the lattice data 
of Ref.[^ for /3 = 5.7, L = 80 (circles) and L = 64 (triangles). 
The standard PT one-loop behaviour is reported as a dotted 
line. 


g 

rr? 

A (GeV) 

m (GeV) 

M (GeV) 

4.0 

0.238 

1.24 

0.60 

0.56 

4.4 

0.274 

1.15 

0.60 

0.56 

4.9 

0.313 

1.05 

0.59 

0.56 


Table I: Masses for = 3 at the energy scale of Ref.Q. 


agreement is fairly good but less remarkable than for the 
gluon, with a deviation that increases as p approaches the 
cutoff. Since it is well known that a fixed-order calcula¬ 
tion needs renormalization group (RG) corrections in the 
UV limit, we do not bother about the large-ener^ pre¬ 
dictions of the model. As largely discussed in Ref.[2|, the 
mass parameter becomes irrelevant in the UV and the op¬ 
timized perturbative expansion is equivalent to the stan¬ 
dard PT in that limit. On the other hand, in the IR the 
optimized perturbative expansion open the way to a first- 
principle perturbative description of the phenomenology 
that can be improved order by order and reach a degree 
of accuracy comparable to the UV limit. In fact, by a 
MOM-Taylor schemefl^ with the vertex renormalization 
constant set to one in Landau gauge, a running coupling 
can be defined by the RG invariant quantity 

a{p^) = ^Aip)F^{p) ( 11 ) 

that has been measured on the lattice in Ref.Q. The 
coupling dose not get too large in the IR and can be eval¬ 
uated in PT, order by order, from first principles, without 
any information from the lattice. In fact, renormalizing 
at the same scale of Figs. 3 and 4, we can evaluate the 
running coupling directly from Ea. (|ll|l and compare the 
result with the lattice data. As shown in Eig.5 the agree¬ 
ment is very good in the IR. 

These findings are encouraging and we expect that 
eventually, including quarks in the formalism would lead 
to a full perturbative description of QCD in the IR. 
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